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THE CYGNUS REGION: 26 AL FROM OB ASSOCIATIONS 
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ABSTRACT 



The COMPTEL map of the 1.809 MeV gamma-ray 
line, which is attributed to the radioactive decay of 
26 Al, shows significant excess emission in the Cygnus 
region. Cygnus is a region of recent star formation 
activity, which is rich in massive, early-type stars 
concentrated in OB associations. We perform pop- 
ulation synthesis studies of this region to model the 
production and distribution of 26 Al, which is pre- 
dominantly synthesized in massive stars and dis- 
tributed through stellar winds and core-collapse su- 
pernovae. 

Our simulations of the Cygnus OB associations de- 
termine the time-dependent production rate of 26 Al. 
We also estimate the associated production of ra- 
dioactive 60 Fe, which is primarily synthesized in core- 
collapse supernovae. Gamma-ray line emission from 
this isotope has yet to be detected. We also track 
the temporal evolution of the mechanical " luminos- 
ity" due to the kinetic energy of stellar winds and 
supernova ejecta, and the ionizing EUV emission 
from early-type stars. Together with a simple one- 
dimensional model of the superbubble driven into the 
ISM by this energy flux one can constrain key param- 
eters of the stellar associations in the Cygnus region. 




Figure 1. Close-up of the Cygnus region from the 
most recent COMPTEL 1.809 MeV maximum en- 
tropy map. The deduced fluxes for different areas 
are given in the insert. 
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1. THE CYGNUS REGION 



Recent COMPTEL maps of the galactic 1.809 MeV 
em ission show enhanced emission in the Cygnus re- 



gioi i. Figure 1 displays a detailed image of this region 
on the sky d erived with the maxim um entropy imag- 
ing method (Pluschke et al., 2000). The flux peaks 
strongly in the area labelled Cygnus West and ex- 
tends over ~ 25°. 

E xtensive comparisons of the 1.809 MeV map to po- 
ten tial tracers ( Knodlscder et al., 1999 ) strongly sup- 



ort t he idea that massive stars (e.g. jMcynet et al 
1997|)) and their sub-sequent core-collapse super- 



novae (e.g. |Woosley fc Weaver (1995| )) are the dom- 
inant sources for interstellar ZD A1 . In addition, the 



studies of Kolb (2000 ), based on detailed modeling of 
binary star populations, suggest that the integrated 
contribution from novae to the galactic 26 Al budget 
is most likely less than 0.1 M Q , i.e., less than 10% of 
the total amount inferred from the 1.809 MeV map. 



A comparis on with other determinations (Starrfield 



et al., 1993) renders the estimation of the nova contri- 
Tbi 



bution to be seriously affected by theoretical uncer- 
tainties. Nevertheless, due to the much longer stellar 
evolution the nova contribution could be neglected in 
our attempt to model OB assoc iations. The contri- 
bution of 26 Al from AGB stars ( Mowlavi fc Meynet, 



It 



2000 ) is expected to be less than ~ 0.4 M Q . 
presently impossible to distinguish between contri- 
butions from hydrostatic burning, injected by stellar 
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winds, and 26 Al from explosive burning taking place 
in the shocks of core-collapse supernovae. Only the 
strong correlation of the observed 1.809 MeV emis- 
sion with tracers of massive stars (most significant 
for free-free emission) provides some clues on the rel- 
ative importance of the two production channels. 
The Cygnus region contains 23 k nown Wolf-Rayet 
stars (van der Hucht ct al., 2000), of which 14 are 
classified as WN stars, 8 as WC and 1 as WO type. 
In addition, recent supernova re mnant (SNR ) cata- 
logues list 19 SNR in this area ( preen, 200"o| ). Fur- 
thermore, 7 OB associations have thus far been as- 
sociated with the Cygnus region. 
We have extended traditional stellar population syn- 
thesis models by including the nucleosynthesis of 
26 Al and 60 Fe in massive stars and supernova ex- 
plosions. This allows us to determine the amount 
of radioactive material as a function of time for a 
given star formation history. In addition to calculat- 
ing the budget of radioactivity in the star forming 
region, population synthesis models also provide the 
energy production rate due to the kinetic energy in 
stellar winds and expanding supernova ejecta, as well 
as the flux of ionising radiation from hot, massive 
stars. It is thus possible, in principle, to reconstruct, 
or at least constrain, the star formation history of the 
Cygnus region by combining data on the dynamics, 
radiation environment, and gamma-ray line fluxes. 



2. THE OB ASSOCIATION MODEL 



Our population synthesis model is based o n the idea 
of continuo us parameter distributions (Pluschkc 
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Figure 2. Evolution in time of the amount of 26 Al 
for various assumptions about the IMF slope and the 
SFR(t) profile. 



wind component, whereas the second peak is due 
to the delayed supernova contribution. For very 
flat IMFs the profiles are very strongly peaked in 
the early phase. Long duration episodes of star 
formation (keeping the total amount of stars pro- 
duced constant) result in flat profiles (essentially a 
situation of steady state) with significantly reduced 
values for the maximum amount of radioactivity. 

Through their enormous mass loss rates, in 
particular during the Wolf-Rayet phase, massive 
stars release a large amount of kinetic energy into 
the surrounding medium. Sub-sequent supernovae 
contribute a large amount of kinetic energy through 
the ejecta moving with initial speeds in excess of 
10000 km/s. By applying the Geneva mass loss 
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description ( Meynet et al., 1997), a semi-empirical 



veloc ity law ( [Howarth fe Prinja, 1989 ; Prinja et al., 



by [PlLlsdikL 1 cL al. (2000 ). Based uii biro assmiipLioii 



and W ^l a ;£ pll ^ d C " Ch ° . , C yS nus re S? on 1990D , and assuming the canonical 10 01 erg per 



of a time-invariant initial mass function (IMF) and a 
mass- invariant star formation rate function, SFR(t), 
we compute models for OB associations of different 
richness, duration of star forming activity, and slope 
of the initial mass function. In all cases the IMF 
is assumed to be a single power-law with mass 



limits Mi, 



8M Q and M u 



120 M Q , and the 

(i+r) 



parameter T determines the slope via $ oc M Z ^ MS 

For the computation of the time-dependent 
ejection of 26 Al and 60 Fe due to stellar winds and 
core-collapse supernovae we use the predicte d yields 
from the detailed nucleosynthesis models of Meynet 
|st al. (1997j ) and |Woosley fe Weaver (1995| ). For 
a given star formation history the production rate 
of 26 Al and 60 Fe follows from stellar life times and 
the time- and mass-dependent yields. Radioactive 
decay then predicts the emerging gamma-ray line 
flux together with the accumulated mass of 26 Al 
and 60 Fe , which is displayed as a function of time 
in figure for the case of 26 Al . The displayed 
profiles are normalised to the total mass converted 
into stars in the given mass range, varying the slope 
of the IMF and SFR(t). As long as the IMF is not 
too flat, the time profiles of populations with short 
duration star formation are characterised by two 
peaks. The first peak is attributed to the stellar 



supernova explosion, we calculated the evolution 
of the mechanical luminosity of an OB association. 
The kinetic energy ejected by an association can 
in principle be observed in form of an expanding 
supershell around the association. Several processes 
coupled to the shell expansion modify the observable 
kinetic energy residing in these shells. To link the 
observed kinematic properties of these shells to the 
underlying energy input rate we extended our popu- 
lation synthesis model by an one-dimensional bubble 
expansion model. This model uses the thin shell 



approximation (Shull & Saken, 1995) and includes 
radiative cooling and mass loading of the bubble 
interior due to partial shell evaporation. We find 
that the kinetic energy derived from observations of 
the shell is typically less than 20% of the injected 
energy. However, this fraction depends strongly on 
the details of the poorly understood evaporation 
process taking place at the interfaces between shell 
material and the hot bubble interior. More efficient 
mass loading of the bubble interior increases the 
density, which leads to more efficient radiative 
cooling. Efficient shell evaporation thus causes 
the expansion to stall sooner. Figure || displays 
the evolution of the relation between the activity 
(decay rate) of the radio isotopes and the kinetic 
energy of the shell. Both quantities can be obtained 
observationally, thus providing a strong constraint 
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Figure 3. Activity (decay rate) of 26 Al & 60 Fe in 
a star forming region and the associated kinetic en- 
ergy of the expanding supershell, driven mostly by the 
stars that are also the producers of these radioactive 
isotopes. 



Fits to tracer maps (Knodlseder et al., 1999) 



have shown that free-free emission from the ionised 
component of the ISM is most strongly correlated 
with the 1.809 MeV map. This motivates the 
inclusion of the emission of photoionizing EUV 
radiation from massive, young stars in our popu- 
lat ion synthesis mod el. Based on the calculations 
of |Vacca et al. (1996| ) we used a simple fit function 
to compute the integrated emission rate. Our 
simplified description overestimates the EUV emis- 
sion by up to 0.3 dex because the effects of stellar 
evolution are not included rigorously. One of the 
main features in the evolution of the ionizing flux 
is its rapid decrease, due to the fact that only the 
most massive stars contribute significantly to this 
radiation. Observations during the early phase of a 
starburst thus provide valuable constraints on the 
stellar population near the upper end of the IMF. 
This is also the regime probed by gamma-ray line 
observations. 



3. CYGNUS OB ASSOCIATIONS 



Based on the OB associations listed in Bochkarev 
[Sitnik (1985 ) we applied our OB association model 
to the gamma-ray observations of the Cygnus region. 
We attempt to constrain model parameters for this 
region by directly using star count data, as well as 
age estimates and shell energy determinations from 
the literature. For the 7-ray intensity model of the 
Cygnus region we assumed a circular shape of the 
emission areas with a central plateau ranging over 
2/3 of the angular radii of the associations and a 
sharp Gaussian drop-off near the border. Fitting of 
thi s OB association model with six independent asso- 



ciaiions to the COMPTEL 1.8 MeV data attributes 
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Figure 4- Comparison of the reconstructed maximum 
entropy image of the association (image) with the 
COMPTEL 1.8 MeV map (contour). 



this model the remaining OB associations contribute 
less than 25% to the total flux. Figure f| shows an 
overlay of the maximum entropy image (contours; 
see Fig. 0) and a maximum entropy deconvolution 
reconstruction of the fitted model (image). 
This result is consistent with the recent work of 



Knodlseder (2000[ ), who found that Cygnus OB2 is 
the richest OB association known to date. B y us- 
ing the 2MASS survey data |Knodlseder (2000| ) con- 
cluded that Cygnus OB2 resembles more closely to 
a young globular cluster, such as the Arches cluster, 
than a typical OB associations. He found about 120 
O star and 2600 OB star members in an area of 2° 
in diame ter. The large number o f OB stars revives 
the idea flComcron fc Torra, 1994|) that Cygnus OB2 



triggered the star forma tion activity in Cygnus OBI, 
OB3, OB7 and OB9. pomeron et al. (1998| ) con- 
cluded from kinematic arguments that Cygnus OB2 
must provide a steady-state mechanical luminosity 
of ~ 4.7 • 10 39 ergs _1 , to get this scenario working. 
This implies that 2800 early type stars should reside 
in Cygnus OB2, which is in ag reement with the ob- 
servation of Knodlseder (2000 ). 
However, if one regards results from an the anal- 
ysis of a c ontinuum radio survey of the Cygnus X 
region by Wcndker et al. (1991), which results in 
a none detection of new supernova remnant inside 
Cygnus OB2, a contradiction arises. G78. 2+2.1 ap- 
pears to be the only SNR in this field, which in turn 
implies that Cygnus OB2 is significantly younger 
than 4.5 Myr. This conclusion i s supported by the 
large number O stars found by Knodlseder (200C ) 



and b y the HR-diagram extracted by Massey et al 



the strongest contribution to Cygnus OB2 followed 
by a 50% smaller contribution from Cygnus OBI. In 



(1995 ), which gives an deduced age of 2.5 ± 1.0 Myr 
for Cygnus OB2. The absence of supernova activity 
inside Cygnus OB2 forces the number of early stars 
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to be four times higher to fulfill the energy criterion 
for propagating star formati on, which is in con tradic- 
tion to the observations of [Knodlscder (2000 ). The 
implications of our model fit (that Cygnus OBI may 
provide a larger amount of 26 Al ) supports this argu- 
ment due to the fact that OBI must be old enough 
to release 26 Al , which is not the case in the propa- 
gation scenario assuming an age of less than 4 Myr 
for Cygnus OB2. Therefore we are inclined to dis- 
regard the idea of propagating star formation by an 
expanding bubble structure around Cygnus OB2. 
However, supernovae appearing inside hot, thin cav- 
ities may not produce visible remnants so the con- 
strain on OB2's age may not be very stringent. If 
ind eed OB2 is older than previously assumed the 
pro posed scenario may work. 



4. UNCERTAINTIES 



Our model is subject to two different types of un- 
certainties. First of all, all of the input components 
of the model have systematic uncertainties in the 
theoretical descriptions used to build the basis of the 
model. For example, in the case of the nucleosynthe- 
sis part, the employed stellar evolution models are 
still hampered by an incomplete understanding of 
internal mixing processes in stars. The description 
of convection has been intensively debated over 
recent years. In addition, uncertainties in the mass 
loss description as well as in the rotation treatment 
strongly affect the yield estim ates . Recent inves- 



tigat ions of stellar evolution (Meynet & Maeder, 
000 ) and sub-seq uent core-collapse super novae 
from rotating stars (Woosley & Heger, 1999) have 
shown stellar rotation to be a key ingredient in 
understanding the coupling of internal mixing and 
mass loss. One of the most severe changes due to 
rotation is a extension of the main sequence phase of 
up to 25% compared to non-rot ating models (using 
same input physics). However, Meynet fc Maodcij 
concl ude that their earlier models (Meyne 



:2000|) 



:t al., 1997) using Schwarzschild convection with 
moderate overshooting together with an enhanced 
mass loss may reasonably describe the behaviour 
of massive stars rotating with a mean velocity and 
avoiding overshooting and mass loss enhancement. 
Nevertheless, the extended main sequence lifetime of 
rotating stars may alter the shape of the predicted 
ISM abundances by shifting and expanding it in 
time. In addition, in some cases the nuclear input 
to the employed reaction networks is not sufficiently 
well known to derive accurate yields. In summary, 
each imperfection of the input models may give rise 
to individual uncer tainties of up to a factor 2. 
In the case of 26 Al Langer & Braun (1998) demon- 
strated that 26 Al synthesis can be strongly affected 
by the presence of a binary companion. Their 
simulations suggest possible yield enhancements by 
factors up to 1000. Clearly the interpretation of the 
observed 1.809 MeV emission in the Cygnus region, 
and the Galaxy as a whole, depends on our un- 
derstanding of this binary star "correction factor". 



More stellar evolution studies are clearly needed in 
this area. We assume that significant uncertainties 
are induced by our lack of a full understanding of 
nucleosynthesis in binary star systems, but it is 
presently not possible to obtain reliable estimates of 
these uncertainties. 

Because stellar mass loss i s still poorly understood 
(Lamers & Cassinelli, 1999), predicting the mechan- 
ical luminosity is uncertain as well. 
Also the determination of the initial mass func- 
tion and the star formation history is subject to 
uncertainties. For our standard case we are using 
the Salpeter-IMF. |Knodlscdcr (2000[) determined a 
steeper IMF slope of T = 1.6 ± 0.1 for Cygn us OB2 
which repre sents a slope as favoured by Kroupa 



et al. (200q ) 



In addition to these systematic uncertainties 
statistical fluctuations due to the finite population 
size of real OB associations may affect the predicted 
light curves. To include this aspect in our model 
we use a Monte Carlo version of our population 
synthesis code. 
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Figure 5. Dependence of spread of the interstellar 
Al mass as function of the richness (number of 
stars more massive than 8 M Q ) of the association. 



Figure a shows the dependence of the simulated 
spread of the amount of 26 Al (measured at the time 
of peak 1.809 MeV emission) on the richness of the 
association. For populations containing more than 
100 stars the statistical uncertainties become less 
than the systematical ones. 



5. SUMMARY & PROSPECTS 



We have added the production of radioactive 26 Al 
and 60 Fe to a population synthesis code in order to 
study gamma-ray line emission from Galactic star 
forming regions. The numerical simulation of a star 
burst (or a more extended period of star forming 
activity) allows us to calculate time-dependent 
gamma-ray line fluxes, which can be compared to 
observations to constrain key parameters, such as 
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the strength, duration, and shape of the SFR(t) 
function. However, the model presented here is 
specifically designed for applications to OB associa- 
tions, i.e., star formation involving a rather limited 
number of stars. We therefore had to investigate 
statistical uncertainties with Monte Carlo methods. 
For OB associations containing fewer than ~ 100 
stars statistical uncertainties dominate. 
Application of our model to the Cygnus region 
demonstrates the use of gamma-ray observations as 
a complementary tool in the study of star formation 
activity in the Milky Way. Subtracting a contri- 
bution of approximately 40% of the total observed 
flux due to isolated sources such as WR stars and 
SNRs we find Cygnus OB2 to be the dominant 
source of interstellar 26 Al in this region contributing 
about 50% of the integral gamma-ray line emission 
attributed to the Cygnus associations. Cygnus 
OBI contributes up to 25% whereas the remaining 
associations altogether add another 25%. Our 
population synthesis disregards the propagating star 
formation hypothesis due to some contradictions 
found in the expected gamma-ray line flux and the 
kinetic energy output. 

Future studies with our population synthesis 
tool will utilize constraints from the observable 
kinetic energy in an expanding supershclls. To 
take advantage of this particular diagnostic tool, 
a careful study of the effects of shell evaporation 
needs to be carried out. We plan to use existing 
HI surveys to obtain data for a larger sample of 
Galactic supershells, to determine properties of 
these shells, and to progress to population synthesis 
simulations of global star formation in the Galaxy. 
Based on the interpretation of existing gamma-ray 
line observations, such as the COMPTEL 26 Al map 
at 1.809 MeV, we will simulate the performance of 
the INTEGRAL mission, especially with focus on 
the gamma-ray lines from 60 Fe that have yet to be 
detected. 
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